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We present the detailed analysis of genuine twist-3 contributions to the exclusive electropro¬ 
duction amplitude of transversely polarized vector mesons. Using the formalism based on 
the QCD factorization in the momentum representation we calculated all the genuine twist-3 
terms and found the total expression of this amplitude at 1/Q level. Generally speaking, these 
terms violate standard factorization owing to the existence of the infrared divergencies in the 
amplitude of hard sub-processes, although the strongest divergencies cancel due to the QCD 
equations of motion. We discuss the possible treatment of surviving divergencies. 

I. INTRODUCTION 

Hard exclusive reactions provide important information for unveiling the composite structure of hadrons. More¬ 
over, a self-consistent description of the hard exclusive reactions is one of main goals of the QCD application. 
The perturbative theory allows to implement the systematical calculations providing that, first, it is possible to 
single out the amplitude of the short-distance sub-process and, second, to prove the infrared finiteness of this 
amplitude. Technically it corresponds to the factorization of large (hard) momenta from small (soft) momenta 
domains. Besides, the factorized pieces of amplitude depend on the dynamics that is typical at a given scale, i.e. 
the amplitude becomes the convolution of hard and soft parts. However, the electroproduction of transversely 
polarized vector mesons 


hadron(pi) -|- "f*{q) p{p) + hadron(p 2 ) , (1) 

provides the well-known example of QCD factorization breaking. Indeed, the factorization is valid only in the 
case of longitudinally polarized p-meson production [1,2], when the factorization assert that the photon-parton 
hard, perturbatively calculable, sub-processes are parted from the nonperturbative matrix elements. 
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Unfortunately, it is not so encouragingly simple for the case of transversely polarized meson production. Its 
description is strongly complicated due to the existence of infrared divergencies in amplitudes, breaking down 
the factorization ( see e.g. [3], [4] and references therein). 

The amplitude of transverse vector meson production corresponds to the contributions suppressed as 1 jQ in 
comparison with the longitudinal vector meson case [1]. At the same time, the recent experiments show that 
the transverse mesons production amplitudes provide a sensible contributions even at moderate virtualities 
[5]. So, to describe these processes one should take into account the l/Q terms. 

In Ref. [2], the analysis of twist-2 amplitudes for hard exclusive electroproduction of mesons in terms of gen¬ 
eralized parton distributions (GPD) was presented. Later, the authors of [3] have discussed the factorization 
problems in the electroproduction of light vector mesons from transversely polarized photons. They have 
taken into account the kinematical twist-3 terms within the Wandzura-Wilczek (WW) approximation. Also, 
the helicity flip amplitude of transversely polarized vector mesons production was considered within the same 
approximation [6]. 

Although the kinematical and dynamical (genuine) higher twists contributions are, generally speaking, in¬ 
dependent, there are notable exceptions in some kinematical regions. In deep inelastic scattering at xb 1 
the kinematical higher twist terms, described by Nachtmann variable, lead to inconsistencies unless the genuine 
higher twists are taken into account [7]. One cannot exclude, that the genuine higher twists may cure, at least 
partially, the problem arising from the treatment of the end-point regions in hard electroproduction. 

Thus, in this paper our goal is to study the role of genuine twist-3 contributions in the factorization theo¬ 
rem breaking. We adhere the approach based on the momentum representation, the basic stages of which 
are expounded in previous papers [8]- [11]. The essence of this approach constitutes of the generalization of 
Ellis-Furmanski-Petronzio factorization scheme [13]. These authors have considered the twist-4 effects in the 
processes of unpolarized deep inelastic scattering. At the same time, our formalism is more closely following the 
approach [14] developed by A.V. Efremov and one of the authors, where the role of twist-3 terms in polarized 
deep inelastic scattering was studied in detail. 

In this paper we have computed the total expression for transverse meson production amplitude comprising the 
quark and gluon leading twist GPDs in nucleons and the genuine twist-3 contributions to the p-meson wave 
function. As a cross-check, we have reproduced the gluon contributions to the transverse meson production 
amplitude in WW approximation obtained in [3]. We also have discussed the possible ways to treat the infrared 
singularities and their partial cancellation. 

The structure of paper is as follows: in Sections II and III we will introduce the kinematics and parameterization 
of matrix elements. Also we will present the total amplitude of quark and gluon distribution diagrams in Sections 
IV and V. We make the conclusion in the last section. 

II. KINEMATICS 

Let us start with the description of our kinematics: the p is momentum of transversely polarized p-meson 
and its polarization vector is e^; the momentum of virtual photon is denoted by q (Q^ = —q^). As the hard 
exclusive production kinematics is similar to the DVGS kinematics, it is natural to use the analogous notations 
(cf. [8]). By making use of initial (pi) and final (p 2 ) nucleon momenta we construct the average momentum P 
and transferred momentum A: 
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( 2 ) 


p ^ P2 +Pl 
2 


A = P2-P1, A‘^ = t. 


It is convenient to write up Sudakov decompositions for all the relevant particles. We choose the light-cone basis 
composed by the physical vectors: P, p. Such a choice is possible because in this paper we assume that the 
initial hadron momentum pi and final hadron momentum p 2 are collinear, i.e. —> 0, and p\ = p\ = t = 0^ 
neglecting all the relevant higher twists contributions arising from the nucleon matrix elements. In addition, we 
neglect squares of meson masses. Therefore the Sudakov decomposition up to kinematical twist-3 terms takes 
the form: 


A = -2^P, e = e • 
p = q — A=p-Pn; 


np + , 

n • n = 


P 


1 4g 

P-P 


where we introduced the normalized vectors n and fi. 


(3) 


III. p-MESON MATRIX ELEMENTS AND THEIR PROPERTIES 

We introduce the parameterizations of the p-meson-to-vacuum matrix elements needed for calculation of am¬ 
plitudes (cf. [16]). Keeping all the terms up to the twist-3 order and using the axial (light-like) gauge 

n • A = 0, (4) 

these matrix elements can be written in terms of the light-cone basis vectors (3): 

(p(p)l^(0)7/x^(2:)|0) = pi{y){e ■ n)pf, + 7’3(2/)e^, {p{p)\ij{0)jf,i aj 'tp{z)\0) = (pj{y)pf,e^, (5) 

(p{p)\2p{0)-f^-ff,tp{z)\0) = iipAiy)epai3selpf}ns, (p(p)|V'(0)757a«* '9j i^{z)\0) = i(PA{y)Pii£pap5elppns, (6) 
(p{p)\-ip{0)-ff,gA^{z2)'ip{zi)\0) = ^{yi,y2)PtJ.eJ, 

(p(p)I^(0)757m5^J(^2)V’(2;i)|0) = iJ{yi,y2)Ppepapselppns, 

(p(p)I^(0)7m* iz2)^izi)\0) = ^iyi,y2)Ppep, 

(p(p)I^(0)757m* Dp iz2)^izi)\0) = iJ{yi,y2)PpepapseaPpns 
^(j/1,2/2) = ^{pl{yi) + •pl{y2^5{yi -y2)P^{yi,y2), 

J{yi,y2) = ^(^PAiyi) + - y2) + J{yi,y2) ( 7 ) 

■« ^ 

here dp = -^{dp — dp) is the standard antisymmetric derivative and = denotes the Fourier transformation with 
measure (zi = Xiu): 

dye~^^^^ for quark correlators, 

dyi (ii /2 for quark — gluon correlators. (8) 

Note that the function pi corresponds to the twist-2; the functions correspond to the WW twist-3, 

functions <I> and J - to the genuine (dynamical) twist-3, while functions p^, pA-, 4*, J contain both. 
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In (5)-(7) the functions (pi, (^3 and ip a parameterizing the two-particle correlators obey the following symmetry 
properties: 


<Pi{y) = ipi{l-y), <p3{y) = <p3{l-y), ipA{y) = -‘PAil - y)- (9) 

At the same time, the symmetry properties of the functions ^ and J parameterizing the three-particle correlators 
are: 


^(2/i,y2) = ^>(1 - y2,i - yi), J{yi,y2) = -J{i-y2A-yi)- (lO) 

The relations (10) represent the particular case of the relations for the functions parameterizing the analogous 
matrix elements in the production of pion pair with the arbitrary angular momentum j [ 10 ], while the case 
under consideration corresponds, obviously, to j = 1. In the case of arbitrary j the two-body functions depend 
on the extra parameter ^ = (p^ —p'^)-n characterizing the skewedness and they have the following properties 

$""(yi,y2;e) =<l>""(l-y2,l-yi;-0, T""(yi,y2;e) =- y2, l - yi;-C- (H) 

For the sake of comparison, we show the relations between our parameterizing functions introduced in (5) and 
functions used in [16]: 

1 (v) 

‘Pliy) ^ fpmp(l>\\{y), Paiy) fpmpg^'"\y), ipA{y) ^--fprUp ^^—. (12) 

Let us remind that the functions psiy) (y^’'^(y)) and pA{y) {g^^Hy)) can be expressed through the pi{y) (<(>|| (y)) 
owing to the WW-type relations ]10] ( ]16]). The WW approximations take the especially simple form in terms 
of functions ] 10 ] 

p±{y) = P3{y)±PA{y)- (13) 

WW relations for p meson represent the particular case of these relations for pion pair with arbitrary angular 
momentum [ 10 ] and take the form: 



X 0 


Let us turn to the integral relations arising from the QCD equations of motion, which may be derived closely 
following [ 8 ] and [ 10 ]: 

1 

J dy(^¥^\x,y) - - 0 psix) + ^Pa{x), 

0 

1 

J dy(^^‘''^\x,y)-J^^\x,y)^=-(x-^^PAix)-^P3{x), (15) 

0 

where symmetric and anti-symmetric functions, whose appearance is crucial for the compatibility with the 
symmetry properties (9, 10), are defined as (/ = $, J): 


^Here we correct the misprints in [10], formula (48). Notations: and correspond to B and D, respectively 
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( 16 ) 


{x, 2 /) = ^ {f{x, y) ± f{y, x)). 

Moreover, the symmetry properties make two equations equivalent to the single one, which may be rewritten 
in the simple form (c.f. [14]): 

1 

J dy 2 T-{yi,y 2 ) = {1 - yi)(p-iyi). (17) 

0 

Here and below we use the obvious similar notations for quark-gluon correlators: 

^±{yi,y2) = J{yi,y2) ±^{yi,y2), :F±{yi,y2) = Jiyi,y2) ±^{yi,y2)- (18) 

As we will see in the next section, the QCD equations of motion play the important role in the cancellation of 
some leading infrared divergencies of the transverse p-meson production amplitudes. 


IV. AMPLITUDE OF QUARK-PHOTON SCATTERING 


At first we dwell on the computation of the quark contributions to the production amplitude As we would 
like to concentrate on the twist-3 effects in p-meson blob in Fig.l, it is sufficient to keep the twist-2 GPD in 
the parameterization of quark matrix elements over nucleon states. Hence, taking into account (2) and (3) the 
parameterization acquires the following form 


{N{p2)\'ip{0)j^^{z)\N{pi)) = H{x)U{p2)yi,U{pi) = - ^'^H{x)Pf, 


(19) 


where, as it was above, = denotes the Fourier transformation, except that the measure now is {z = Ah) 

^20) 

Further, using the parameterization (5)-(7), (19) and calculating the traces, the amplitude given by the simplest 
Feynman diagrams, Figs. 2a and 2b, reads (cf. [2]): 

„T rl 




(q),iT^PT _ o /i _ p2 f/i 

°V ^ ? Ar /n2 


Nc J_ 


dxH (x) 


1 


1 


X — ^ + ie X + ^ — ie 


S. 


( 9 ) 


where 
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(9) 


Jo y 


( 21 ) 


( 22 ) 


Note that the double poles in x cancel due to the use of the gluon propagator in the axial gauge. 

It is convenient to organize genuine twist 3 diagrams according to the insertions of extra gluon to the diagrams 
of Fig.2. Summing the diagrams of Figs. 3 and 4, we finally obtain the expression for the quark distribution 
amplitude that includes all the WW and genuine twist-3 contributions: 


^(9b 77 -Pt = 8y/rVe J j X J) 


(9) 

2 


( 23 ) 


^Recently we computed this amplitude using another parameterization [12] 
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where 


Til ~ ^ J dxH{x) 


1 1 
+ 


, Ti 2 = J dxH{x) 


X — ^ + ie X + ^ — ie 


and 


= J dyi dy2 1 


(x + ^ — (x — ^ + J ' 

4CF^_(yi,2/2) , CAd^-{yi,y2) , 2CF^_(yi, ^ 2 ) - ^ 2 )' 


(1- 2 / 1)2 (l_y2)(l-2/l) 


yiil + yi - 2 / 2 ) 


(24) 


(25) 


= [ dyi dy2 

Jo 


CFiF-iyi,y 2 ) 2 Cf.F_( 2 /i, 2 / 2 ) - Ca-^-( 2 / 1 , 2 / 2 ) 
(1 - 2 / 2 )(l - 2 / 1 ) (1 + 2 / 1 - 2 / 2 )(l- 2 / 2 ) 


(26) 


Note that the entire result is expressed in terms of ” — ” combinations of twist-3 functions (which may be 
transformed to ” + ” ones by making use of symmetry properties (10)). Such combinations of axial and vector 
twist-3 matrix elements are known to appear also in inclusive case [14,17,18]. 

The Til- and -structure integrals possess the poles of second order. This fact leads to the vulnerability 
of factorization theorem. However, one can see that the first term in (25) may be reduced to the two-particle 
function (/?_( 2 /i) by making use of the QCD equations of motion (17). Asa result, the power of infrared divergence 
is reduced due to the factor 1 — 2/1 in its r.h.s. Moreover, in the case of the gauge different from the axial one 
(Feynman gauge, in particular) the emerging double pole in a; in (21) is cancelled by the corresponding additional 
contribution to (23) (together with the dependence on the gauge fixing parameter) provided the equation of 
motion are taken into account. This is quite natural: while the longitudinal amplitude is gauge invariant by 
itself, the gauge invariant set of diagrams for the transverse amplitude contains also quark gluon diagrams, 
related to the quark ones by the equations of motion. 

Still, the remaining single poles would lead to infrared divergencies unless functions $, $, J, J vanish at //i ^ 1 
or 2 /i ^ 0 at least as the first power of (1 — 2/i) or yi. This behaviour is quite reasonable for genuine twist 
term. However, as the corresponding integral (with a colour factor Ca) becomes finite providing they have 
the corresponding behaviour at the edge points, it cannot be used to cancel the infrared divergence of WW 
contribution in the term with the factor Cp. 

To be more specihc, let us calculate this WW contribution explicitly. Within the WW approximation where all 
the quark-gluon parameterizing functions +( 2 / 1 , 2 / 2 ) and ^{yi,y 2 ) (see, (7)) are equal to zero, the integrals at 
the structures Tii and 0.2 in the sum of (21) and (23) take the following forms: 

One can see from (14) that the simple pole in of (27) is cancelled by the zero boundary condition for the 
function ip-{y) : V5-(0) = 0. This is quite similar to the cancellation providing the finiteness of the two pions 
production amplitude in the collision of real and virtual photons [10]. This makes finite the coefficient of the 
potentially dangerous (due to the double pole in x) integral TLi. At the same time, there is no such effect for 
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V. AMPLITUDE OF GLUON-PHOTON SCATTERING 


In this section we turn to the calculation of gluon distribution amplitude. Since we will deal with the two-gluon 
GPD, it is natural to choose the gauge-fixing condition for the gluons in the form: 

n • A = 0, (28) 

where n-vector (3) corresponds to the ’’minus” component of the light-cone basis where ’’plus” component is 
provided by nucleon (average) momentum. In this case, the gluon fields can be expressed through the 
gauge-invariant field-strength tensor G^jy. Consequently, one can parameterize the nucleon matrix elements of 
two gluon fields in terms of the gauge-invariant gluon distribution: 

(JV(b)|.4:(0K(J)I«(m)) £ - P.n, - ^ ^ ^ ^ . (29) 

Note that the gauge condition we use here is different from the one used for parameterization of vector meson 
matrix elements (see, (5)-(7)); so the question arises, is it possible to use these expressions here. Therefore, let 
us discuss this problem in more detail. 

In principle, owing to the Lorentz and gauge invariance, our physical amplitudes should be independent on the 
explicit choice form of n and fi. In the general case, the vectors n and fi can be chosen in an arbitrary way [10]. 
Now we would like to make the comparative analysis of two gauges: (4) and (28), used for p-meson matrix 
elements. In the case of gauge (4) that was also used for the quark-gluon scattering calculation, the vector n in 
(5)-(7) play the double role. Namely, it fixes the gauge and defines the longitudinal momentum fraction x = kn 
carried by the active quark with the momentum k in p-meson. At the same when the gauge (28) is adopted, 
these two roles are distributed between the vectors n and n. Namely, n defines the longitudinal momentum 
fraction whereas fi defines the gauge-fixing condition. Of course, in this case the parameterization of relevant 
p-meson matrix elements should, generally speaking, look more complicated due to the added n-terms. However, 
for ’’physical” n,n (3), one may explicitly check that the parameterization (6)-(7) is still self-consistent. This 
is actually due to the twist-3 approximation, while there is no hope to describe twist-4 terms in such a simple 
way. 

In the gauge (4) only the transverse (physical) gluons exist. In contrast, in the gauge (28) the transverse gluon 
fields are constructed as 


Aj = Ap- ppTi ■ A, (30) 

where the second term in (30) corresponds to the twist-2 contribution. As a result, the use of gauge (28) leads 
to the appearance of twist-2 functions in the parameterization of three-particles p-meson matrix elements, for 
instance 

{p{p)\'>Pi0hpgAp{z2)ipizi)\0) = ^o{yi,y2)ie ■ n)ppPp + ^{yi,y2)ppe^■ (31) 

The new twist-2 terms, proportional to n-A, should be absorbed to the standard P-ordered exponent (the gauge 
link) of the matrix elements of leading twist-2 quark operators [9]. 

Thus, with the help of parameterization (5)-(7) and (29)), we sum the amplitudes given by simplest diagrams 
of Figs. 2c and 2d, with the amplitudes corresponding to diagrams of Figs. 5 and 6. The last sort of amplitudes 
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includes both the kinematical and dynamical (genuine) twist-3 in p-meson-to-vacuum matrix elements. We 
derive the following expressions: the simplest diagrams contributions reads 






1 V 2 


—.fJ dxGix) 


(x + ^ — ie)^ (x — ^ 






dxG{x) 


X — ^ + ie X + ^ — ie 


S. 


(9) 


where 


^>”=1 f ()*>+(:/) 


= 


r dy + 

/o y 


and the genuine twist-3 contributions take the form 


a: 


(9). 7 t^Pt _ 
2,(1 ' 


CfN, 

(iV2 - 1)2 Q 


+ g2xT', 


r(g) 


where 


g.=£dxGw[ (^ + )_.^), + 


and 


n~i9) _ 

-i-l — 


dyi dy2- 


{x - ^ + *e)^J ’ 
CAA-{yi,y2) 


& = 


e7- 


dxG{x) 


1 


1 


= / dyi (i7/2 
Jo 


GpA^ {yi , 2/2) + C'a-^- ( 2/1 , 2/2) 


(32) 


(33) 


(34) 


(35) 


(36) 


'(1-2/2)(1-2/i)’ ^ (1 _ y2)(2/i - 1) 

Note that the integral is essentially non-Abelian and contains only genuine twist-3 contribution. Con¬ 
sequently, the assumption on the linear decrease of T-{yi,y 2 ) leads to its finiteness. At the same time, the 
integral X^'^ becomes divergent due to WW term with the coefficient Gf ■ 

To study these effects and to reproduce the results obtained in [3] we address WW approximation. As in the 
quark case, this approximation implies that all the quark-gluon parameterizing functions J(2/i, 2 / 2 ) and 4*( 2 / 1 , 2 / 2 ) 
(see, (7)) are equal to zero. Using the results of [10] after some computations we derive the following expression: 


where 


^ § I®' }■ 


(37) 




(38) 


If we take into account (12) then the integrals (38) will be rewritten in the forms which completely coincide 
with [3] 


= 




lo y 


2y{l-y)J 

Jo y \ y 22 /( 1 - 2 j 


-y) ' 


(39) 


'^The definition of function >1>|| can be found in [16]. 



Similar to the quark distribution case, the integrals (38) (and, consequently, (39)) are not free from the infrared 
divergencies owing to non-zero boundary values [3,16]. 


VI. DISCUSSION AND CONCLUSIONS 

Let us stress that we calculated the full set of genuine twist 3 diagrams of transverse vector meson electropro¬ 
duction. We also present, for the first time, the WW contribution proportional to quark GPD and reproduce 
the WW contribution due to the gluon GPD calculated earlier by the different method. As a result, we observed 
the cancellation of the few types of the leading infrared divergencies: 

i) The double poles in x of the quark contribution (Fig 2a,b) in the non-axial gauge cancel with the contri¬ 
butions of quark-gluon diagrams (Fig 3,4), provided the QCD equations of motion are taken into account. 

ii) Making use of the equation of motion allows also to eliminate the double poles in y of the same quark-gluon 
diagrams, surviving in the axial gauge, which are reduced to the single poles. 

iii) The single pole in y of the WW contribution proportional to the integral Tii of quark GPD cancels between 
the vector and axial (Pa) distributions. 

Let us start the analysis of potential surviving divergencies from the double poles in x. Gonsider first the 
imaginary parts of TYi, Qi (providing in the gluon case the dominant contribution at large energies), which are 
easily calculated and is represented by the derivatives -^H{x, ^)|a;={, ^G(a;, C)U={ in the transition point x = ^ 
from DGLAP to ERBL region. One may worry, that for the quark case this derivative is not continuous just at 
that transition point due to the existence of so-called Polyakov-Weiss term in the two-component model of GPD 
[21], which is non-zero only for \x\ < However, this term does not give any contribution to the imaginary part 
of amplitude(like any meson exchange term), and should not therefore be taken into account when derivative 
is calculated, so that it should be understood as ■^H{x,^)\x=^+e- Another way of PW term treatment [22] is 
the consideration of its most general form, occupying the whole region \x\ < 1, and being therefore smooth at 
X = ±^. It may be reduced to the original PW form by the sort of ’’gauge transformation” [22], generating the 
irregularities both in PW term and another component of GPD, so that they are cancelled in their sum. The 
possible discontinuity of derivative in the gluon case may be treated by adding the small mass to the quark 
propagators [17], resulting in the similar expression -^H{x,£P)\x=^+e- Anyway, the imaginary parts of 
are well defined, while the real parts may be restored by making use of the dispersion relations [23]. 

Let us now discuss the surviving divergencies arising from the integration over y. The truly non-Abelian higher 
twist contributions {Ca terms) may be assumed finite, as it is sufficient to have the functions ^{yi,y 2 ), J{yi,y 2 ) 
which are going to zero when yi, 2 )j/i ,2 0. The counterpart of this assumption is, however, the impossibility 

to cancel the divergencies of WW terms, which contain another colour factor Cp- 

The only real danger is therefore coming from the WW terms, which have the non-zero values at of (p^, (pA at 
y = 0,1. Even if they are assumed to be zero at some reference point Qq, this property would be (rather slowly) 
broken by the QCD evolution. 

A visible way to regularize these singularities is to keep the transverse momentum kf in loop integrations, 
which is in fact required for the quantitative description of the longitudinal amplitude as well [20]. At large 
the behavior of structure integral is determined by the region where, for instance, I — y is of the same order 
as {k^)/Q‘^. In order to estimate the relevant contribution, staying in the collinear approximation, one may 
implement the corresponding infrared cutoff in the integration over y. As a result, the transverse amplitude 
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is logarithmically enhanced, so that this non-factorizable contribution defines the asymptotic behaviour of the 
ratio of transverse and longitudinal amplitudes: 




niplnQ 


(40) 


where power suppression comes from the standard kinematical enhancement of longitudinal polarization. 

To make the more quantitative estimates one should take into account the factorizable contribution (involving 
the parameterization of genuine twist matrix elements) and include the fer-dependent distributions to non- 
factorizable one. The consistent account for all powers of kx would be equivalent to the summation of all 
kinematical high twist terms. 

To summarize, in this paper we have computed both gluon and quark contributions to the transversely polarized 
p-meson electroproduction up to genuine twist-3 accuracy. We observed a number of interesting cancellations 
of leading infrared divergencies due to the gauge invariance manifested through QCD equations of motions and 
suggested the possible ways of treatment of the surviving infrared contributions. 
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FIG. 1. General structure of factorized electroproduction amplitude 







FIG. 3. Genuine twist-3 diagrams with quark GPD: insertions to the diagrams of Fig. 2a 



FIG. 4. Genuine twist-3 diagrams with quark GPD: insertions to the diagrams of Fig 2b 












